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ficients E^uiE, B) and /i M „(i?, B) is examined. The regime of the eikonal approximation is considered. 
The Hadamard method of field disturbances is used and the dispersion relations are obtained by 
solving the Fresnel equation. Some applications of the formalism are presented. Particularly, bire- 
fringence phenomena induced by applied external fields are derived and discussed. It is shown that 
magnetoelectric birefringence effect can occur even without the presence of Kerr and Cotton-Mouton 
effects, provided the physical system satisfies certain conditions. 
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I. INTRODUCTION 

Inside material media the Maxwell equations must be 
supplemented with constitutive relations between the ap- 
plied external fields and its induced excitations. Such 
relations are generally nonlinear (linear constitutive re- 
lations have also been considered in the literature 
and depend on the physical properties of each considered 
medium under the action of external fields. A remarkable 
consequence coming from the nonlinearity of the field 
equations is the artificially induced birefringence phe- 
nomenon: waves with different polarization propagate 
with different velocities due to the presence of applied 
external fields @, Q. In fact, the effect of applying an 
external field is just to induce an artificial optical axes, 
which in general disappears as soon as the external field 
is turned off. Hence, the initially isotropic medium be- 
comes anisotropic under the action of external fields. In 
the context of crystalline systems the optical axes can 
be present in the medium and birefringence occurs nat- 
urally [in. In this case, by applying an external field 
both effects (natural and induced birefringences) inter- 
fere, leading to a nontrivial mechanism of birefringence 
control by means of external electromagnetic fields Q . 

In the regime of intense electromagnetic fields the 
Maxwell linear theory is no longer applicable and new 
effects, including birefringence, emerge from the quan- 
tum regime \j\ . The analysis of light propagation shows 
that there is a non-null probability of photon splitting 
under a strong external electromagnetic field In- 
vestigations on light propagation in the context of non- 
linear Lagrangian for electrodynamics can be found in 
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fields has been long ago reported in the literature [la . 
EE EE mm m. Distinct proposals for its theoreti- 
cal description were presented in [2l|, HH [HI, H3 > where 
some new aspects of this phenomenon were predicted. 
Among them, the existence of the so called magnetoelec- 
tric birefringence was discussed, which is a kind of bire- 
fringence linear in the product of the electric and mag- 
netic fields. The experimental observation of the magne- 
toelectric birefringence was recently reported in the lit- 
erature dHHi]. 

Birefringence is nowadays widely used in the technol- 
ogy of optical devices as well as a technique for investi- 
gating properties of several physical systems, including 
astrophysical phenomena [271 ] . 

In this work, birefringence phenomena are investigated 
in the context of homogeneous dielectric media at rest 
with the dielectric coefficients £ M „ = e fl l/ (E,B) and 
/i^u = ^^(E, B) in the limit of geometrical optics. The 
analysis is restricted to local electrodynamics, where dis- 
persive effects are neglected. Only monochromatic waves 
are considered, thus avoiding ambiguities with the veloc- 
ity of the wave. 

In Section [IT] the eigenvalue equation associated with 
the light propagation in general non-dispersive mate- 
rial media is presented. In Section IIIII the dispersion 
relations are obtained for isotropic media described by 
= ^u{E,B) and p^ v = ^ V {E,B). Applications for 
birefringence phenomena are investigated in Section lIVl 
where some particular cases are obtained from the for- 
malism. The consequences of the anisotropy present in 
the propagation of the extraordinary ray are discussed. 
Some final remarks are presented in the conclusion sec- 
tion. 

A Minkowskian spacetime employing a Cartesian co- 
ordinate system is used throughout this work. The back- 
ground metric is denoted by ?7 M „ = diag(+l, —1,-1,-1). 
All quantities are refereed as measured by the geodetic 
observer V M = cSq , where 8^ denotes the Kronecker ten- 
sor. For any quantity = (0,A) we define its mod- 
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ulus as X = {—X^X^) 1 / 2 and the associated unit vec- 
tor as X = X/X. In the same way we use the nota- 
tion X^ = X^/X. For any two quantities X^ and 
we denote its scalar product X^Y^ as (XY). As an ex- 
ample, the electric field is represented by — (0, E), 
whose modulus is E = (—E a E a ) 1 / 2 . Similarly, the mag- 
netic field is £? M = (0, B). Its scalar product is given by 
(EB) = E^B^ = -E-B. The units are such that c = 1. 



II. EIGENVALUE EQUATION 



The electrodynamics in a continuum medium at rest is 
completely determined by the Maxwell equations [28| 



V' l D a ^+r 1 a ^%H StP = 0, 



0, 



(1) 

(2) 

together with the constitutive relations D a — s a pE& 
and H a = fi a pB@ . The coefficients e a = E a (E^, B») 
and n a fi = fi a p{E^ , B^) represent the dielectric tensors. 
They are usually denoted as permittivity and permeabil- 
ity tensors, respectively. All the information about the 
dielectric properties of the medium can be encompassed 
in it. 

The propagation of the electromagnetic waves in the 
eikonal approximation of electrodynamics [i| is obtained 
by making use of the method of field disturbances [29l[30l |. 
With the notation introduced before Q we set [i? M ^] s = 
e^K v and [-B'VJj, = ^Kv, where e M and are related to 
the derivatives of the electric and magnetic fields on the 
surface of discontinuity E, and they are associated with 
the polarization of the propagating waves. The quantity 
K\ = dYi/dx x is the wave 4- vector normal to E. Apply- 
ing these boundary conditions to the field equations (fl]) 
and ([2]), we obtain the eigenvalue equation 



Z a T e T = 0, 



where the Fresnel tensor Z a T is given by 



(3) 



1 f)pOi 

C a -I % y"T Ffiv a 



1 dfi \ p A 1 , 2 

Ti?7^ sB VryqpA 2 yl ^ r 

uj oE T 



-q z H x *I a T + H/q q a - H x a q*q s h a T ) , (4) 



with uj = K a V a representing the angular frequency of the 
electromagnetic wave. Additionally, the following defini- 



tions are introduced: 

q a = 
q 2 = 

C a T = 
H a T = 

1 T 



S a T -V a V T , 

h a T K T = K a -ujV 01 , 

-q a q a = lu 2 - K 2 , 

£ T + d& ' 



fl 



dB T 

q a qr 

q 2 



B p , 



(5) 
(6) 
(7) 

(8) 
(9) 
(10) 



The general solution for the wave propagation can be 
derived from the eigen- value problem stated by the gener- 
alized Fresnel equation [Eq. ([3|], and is formally given by 
det \Z a p\ = 0. There are different ways to solve this prob- 
lem. In the next section an expansion of the polarization 
vector in a suitable basis of linearly independent vec- 
tors will be considered in order to obtain solutions from 
the eigen- value equation. 



III. THE DISPERSION RELATIONS 

We shall examine the wave propagation issue in ma- 
terial media with dielectric coefficients depending on the 
intensity of the external fields. For this case we set: 



e a T = e(E 7 B)h a T , 
ii a T = [r l (E,B)h° 



(11) 

(12) 



Since these coefficients do not depend on the directions 
of the external fields, the results are applicable to the 
study of wave propagation in isotropic liquid media. The 
application of the results to crystalline structures can be 
done by considering the explicit dependence of E^ or B^ 
vector fields on the dielectric coefficients. These cases 
were partially described in 0, @| ■ 
From Eqs. (fTTjl and (fl2|) . we obtain: 



de c 



dE T 
dE T 



= -e'h a p E T , 
= +^h a p E T , 



de c 



dB T 
dB T 



= -eh a f3 B T , (13) 



M 2 



= +^h a p B T ,(U) 



where we have defined X' = (l/E)dX/dE and X = 
(l/B)dX/dB, for any quantity X. 

By introducing the last results in Z a p, from Eq. Q, 
it yields: 



Z a T = 



flUJ Z ■■- - 

[ fx 2 



2 2 
jl UJ 



1\ a , M 

- - q qr + -2"^ 



h a T - e'E a E T 



B a B T 



+ 



/' 



fJ?UJ 2 



(qB)B a q T 



fl 



fJL 2 Ul 2 



(qB)q a B T 



-rf^ T E a 



UlfL 



y iSa E T 



qpV^Bs (15) 
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where the following additional notation was introduced: 

(XIY) = X^ V Y V = (XY) + {qX) ^ Y) , (16) 

for any quantities X and Y. 

In the present form, it is clear that Z a T is a 3- 
dimensional object. Particularly we note the Z a T V T = 0. 
In order to find solutions of the eigen- value equation, Eq. 



© with Z a T given by Eq. (|15p . we shall consider the ex- 
pansion of the polarization vector e T in a convenient basis 
of the 3-dimensional space as 

e T = aE T + bB T + cq T . (17) 

In this case, from Eqs.© and (fTS")) . we obtain 



-s'(qE) - 



m 

UJ/I 2 

B a + 



{BIB) + e'E 2 



-V 



an'E 2 



(iq (EB) 

2 2 



qpV^BsEr 

KqB){qE) , fl^E 2 



-e'(EB) 



an'(EB) 



1 ffiB 2 1 



1 ffiB 2 1 



0^ \ /i z 

»B 2 {qB) 

U> 2 [1 2 



-,) {qE) 

Ul/J 2 



2 2 2 

fi(qB)(EB) jfi'E 2 



UJ/l 2 



U> 2 fJ? 



7/Agg) 

LUfM 2 



0. 



(18) 



where a, (3 and 7 are given by 

q(BIB) 



P = 



V / (EIE)(BIB) - (EIB) 2 ' 

-g(£/fl) 
y / {EIE)(BIB) - (EIB) 2 ' 
(gE)(BIB)-(gB)(EIB) 
q v / (EIE)(BIB) - {EIB) 2 



(19) 
(20) 
(21) 



Since {E a , B a , q a } are taken to be linearly independent 
vectors, we obtain the solution to the above system by 
setting each of its coefficients to zero. With a conve- 
nient notation, it yields the following system of algebraic 
equations: 



a{Ai}-6{Bi}-c{Ci} = 0, 
a{A 2 } + b{B 2 }-c{C 2 } = 0, 
a{A 3 }~b{B 3 } + c{C 3 } = 0, 



(22) 
(23) 
(24) 



where 

A! 

A 2 
A 3 



£ + \ - 

i 



UJ/X J flUl 



1 + - (BIB) 



[qVBE] , 



fiq 2 (EIB) (5n'E 2 



U> 2 [X 2 



UJfI z 



-^[B 2 (qE) + (qB)(EB)} 



flW 2 



-ip!E 2 



(25) 
(26) 

(27) 
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(QE) 
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(28) 
(29) 
(30) 

(31) 
(32) 
(33) 



= V aMS qoVpB^Es = q- (E x B) 



with [qVBE] — ,/ Ha y f3 ±J 1 

The above set of algebraic equations [Eqs 
can be solved immediately, resulting in: 



(A 3 B 2 + A 2 B 3 ) d + {A 3 Bx 
+ {A 1 B 2 + A 2 B 1 )C 3 =0. 



A!B 3 )C 2 



(34) 



This is the general equation governing the phenomenon 
of electromagnetic wave propagation inside material me- 
dia described by the dielectric coefficients given by Eqs. 
(fTTj) and (fT2")) . Such equation is usually called as the 
dispersion relation. In the next section some electro- 
magnetic-optic effects will be derived from it and known 
results concerning birefringence phenomena will be re- 
covered and discussed. 

For several physical configurations, the dispersion rela- 
tions can be presented in the suggestive form g± K^K U = 



4 



0. The symmetric tensors g^ represents the optic met- 
rics and, generally, present a solution for each possi- 
ble polarization mode inside the medium. The integral 
curves of the vector AT M are geodesies in the associated 
effective geometry [3j. 

Before closing this section some comments on the ap- 
plicability of the vector basis introduced in Eq. (TTT)) are 
in order. In the next section we shall apply the results 
obtained here to study some limiting cases in which two 
of the vectors used in the basis are parallel. In fact, even 
for those cases the dispersion relation stated by Eq. (|3~4"|) 
holds. It can be understood as follows. Let us consider, 
for instance, an specific configuration where the angle e 
between the vectors E and B is small. Thus, the products 
between these vectors can be presented as £M? ~ 1 — e 2 /2 
and \\E x B\\ ~ e. For £ ^ the set of vectors introduced 
in Eq. (|17[) will still be a set of basis vectors. Now, if 
we take e to be sufficiently small, such that its contri- 
bution in the dispersion relation is not measurable, the 
obtained results must be the same as if E and B would 
be parallel vectors. In this way, the results derived in 
this section can be naturally extended to deal with these 
limiting cases. 

In fact, the dispersion relation obtained in this section 
can also be obtained without taking in consideration an 
specific vector basis. Using the Cayley-Hamilton method 
(see, for instance, Ref. [3l[ ), the eigenvalue equation 
can be solved by means of (Zi) 3 — ZZ\Zi + 2Z$ = 0, 
where Zi (i = 1, 2, 3) represents the traces of tensor Z^ v . 
For simple cases, as occur in the derivation of the Kerr 
effect, the use of the Cayley-Hamilton method is conve- 
nient. Nevertheless, for more elaborate situations, as it 
occurs in the magnetoelectric birefringences, the use of 
this method leads to laborious calculations. As simple 
application of this method the Kerr effect is derived in 
the appendix. 



IV. APPLICATIONS TO BIREFRINGENCE 
PHENOMENA 



tions to the existence of the Kerr and the Cotton-Mouton 
birefringence effects in material media, as described by 
the present formalism. 

For this case, it follows from Eq. (|34j) together with 
definitions stated in Eqs. (|25|) - ([33]) . that: 



Aiw 4 + A 2 v 2 + A 3 = 0, 



(35) 



where we have defined the phase velocity v 2 = w 2 /q 2 , 
and 



Ai= e 2 (s + s'E 2 ), 



(36) 



A, 



A;, 



[B 2 {qEf + {EBf + 2(EB){qE)(qB)] 



e 2 fi(BIB) 2e 2 + ee'E 2 + e e'(qE) 2 



e'fi{BIB)(qE) 2 + e[i{BIB) e + s'(qE) 2 



(37) 



.(38) 



Finally, we obtain the following quadratic equation for 
the phase velocity: 



-A 2 ± y/A 2 - 4AiA 3 
2Ai 



(39) 



In general, we can obtain two solutions from the above 
equation, which correspond to the two possible polariza- 
tion modes propagating in the medium. 



1. Kerr birefringence 

By considering /i = li c — constant in the Eq. (|39|) we 
obtain the following solutions: 



I 



(40) 



and 



In this section some specific examples will be derived 
from the general dispersion relation presented in Eq. 
([34]) . First we will recover some well know cases, as the 
Kerr electro-optical and the Cotton-Mouton magneto- 
optical effects, in which the birefringence phenomena ap- 
pear as a consequence of the applied external fields. After 
we will describe the recently measured [III H|| magneto- 
electric birefringence. In particular, the anisotropic be- 
havior of the extraordinary ray with respect to the direc- 
tion of propagation will be carefully examined. 



1 



^ c {e + e'E 2 ) 



(41) 



where the index o stands for the ordinary ray (o-ray), 
which propagates isotropically, and the index e stands 
for the extraordinary ray (e-ray), which depends on 
the direction of wave propagation. Particularly, if the 
propagation occurs parallel to the external electric field, 
(qE) = q ■ E = 1, the velocities of both rays coincide. 
The difference between them achieve its maximum value 
when (qE) = q ■ E = 0. In this case Eq. (|4"Tj) results in: 



A. Electro and magneto-optical effects 



1 



J e± 



ude + e'E 2 )' 



(42) 



Let us now consider the particular case where e = e(E) 
and \x = li(B). These functional dependence of the di- 
electric coefficients are, respectively, the necessary condi- 



By considering the expansion of the permittivity as 

e = e c + e 1 E 2 1 (43) 
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with e c and e% constants, we obtain 

o o 1 



'II 



fi c e c {l + e 1 E 2 /e c y 
1 



(44) 



(45) 



e± fi c e c {l + 3eiE 2 /e c y 

It is assumed that E\E 2 corresponds to a small correc- 
tion compared to the background permittivity e c . The 
corresponding refraction indexes (n = 1/v) are given by: 

n,| = ±- ~ (l + %^ , (46) 



2e c 



— ~ V£cA« c ( 1 



ie 1 E 2 
2e c 



(47) 



Finally the observable quantity defined as the maximum 
difference between the refraction indexes in the system is 
given, up to second order corrections, by 

, e 1 E 2 , , 

n±_ - n\\ = y/EcVc • (48) 

This result is known as the Kerr electro-optic effect. For 
this case, the presence of an external magnetic field does 
not produce any change in the results. 



2. Cotton-Mouton birefringence 

Now, taking e = e c = constant in Eq. ([55]) . we obtain 
the following solutions: 

v 2 = — 

£c.H 

and 



(49) 



1 

£cH 



l+ l -(BIB) 
A* 



(50) 



As before, the velocity of the e-ray depends on the di- 
rection of propagation. In the particular case where the 
propagation occurs parallel or anti-parallel to the exter- 
nal magnetic field (q ■ B — ±1), its velocity reduces to 
the velocity of the o-ray. The difference between the ve- 
locities of the ordinary and extraordinary rays achieve its 
maximum value when q ■ B = 0. In this case, from Eq. 
USD}: 



J e± 



1-^B 2 



1 

e c fi \ [i 

By considering the expansion of the permeability as 

/! = [lc+ S2 B 2 , 
with [i c and e 2 constants, we obtain 

1 (, Z2B 2 

1 (, ie 2 B 2 



1 



/'<:■ 



(51) 

(52) 

(53) 
(54) 



By considering e^B 2 « fi c , the corresponding refraction 
indexes are given by 



m ~ y/e c n c 1 + 



62 B 2 



n± ~ y/e c n c 1 + 



2fi, 
3s 2 B 2 
2 He 



(55) 
(56) 



Now the difference between these refraction indexes is 
given by 



nj_ — n\\ 



e 2 B 2 

He 



(57) 



This result is known as the Cotton-Mouton magneto- 
optic effect. The presence of an external electric field 
is not important for this case. 

These two effects, Kerr and Cotton-Mouton birefrin- 
gences, are really symmetric effects. This symmetry is 
related to the fact that the Maxwell equations in the ab- 
sence of sources are symmetric with respect to the duality 
rotation. 



B. Magnetoelectric birefringence 

As another application of the present formalism, let 
us now examine some features of the so called magneto- 
electric birefringence. For simplicity, we shall assume in 
this section the particular case where e = e(E, B) and 
u = He = constant. 

In this case, from Eq. (|3~4"1) we obtain the following 
solutions for the phase velocities of both the ordinary 
and the extraordinary rays: 



±- 



1 



(58) 



and 



e [qVBE] 



± 



l e 2 [qVBEY 



-e'(qE) 2 



2(e+e'E 2 ) V A(e+e'E 2 ) 2 ^(e+e'S 2 ) 



.(59) 



As one can see, the equation for the velocity of the e- 
ray contains two possible solutions, which are not obvi- 
ously symmetric with respect to the direction of propa- 
gation. Together with the o-ray we generally have three 
distinct solutions for the phase velocities. In fact the two 
solutions associated with the e-ray indicate that it can 
present distinct velocities in opposite directions. This 
aspect will be further addressed in the analysis of some 
particular cases. 

Let us consider the following expansion for the permit- 
tivity coefficient 



e = e c 



eiE 2 + s 2 B 2 + e 3 (E ■ B) 



(60) 



For practical cases, the terms in £, (i=l,2,3) represent 
just small corrections to the main term e c - the back- 
ground permittivity By using this expansion we obtain 
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the following results for the velocities of the above rays, 
up to second order terms in e.j, 



:±- 



1 



1-^E 2 - 

2e n 



£2 
2e c 



B 



£3 
2e r 



EB(E ■ B) 



and 



(61) 



e 2 



EB 



±- 



£3 

2£ C 
1 4 



E 2 (E ■ B) 
(q-Ef 



v / £cM^ 
-^LEB(E-B) [{q-Ef 



q ■ (E x B) 

- -E 2 

2 



£2 
2e r 



(62) 



Now, two particular cases of current experimental in- 
terest are analyzed. First we shall derive an example 
of the recently measured magnetoelectric linear birefrin- 
gence 26]. After, the so called Jones effect will be pre- 
sented. The latter is also a kind of magnetoelectric bire- 
fringence which was predicted long ago [21( and was re- 
cently measured [25j . 



1. Some specific cases 

Let us consider the case where the direction of the wave 
vector q is perpendicular to the plane which contains the 
electric E and magnetic B fields. Let us work in a Carte- 
sian coordinate system and set E ■ B = cos 9 and q = z 
such that q = (E x B)/\\E x Up to second order 
terms, the velocity of the e-ray, given by Eq. (|62| . is 



P±Q, 



where 



P = — EB + — E 2 cos 9) sin 9, 

' z 2e c 



(63) 



(64) 



Q = 



lJ e ±E 2 -^B 2 - £ -^EBcos9\ .(65) 



V £ cMc V 2 £ c 2e c 



Since e, refers to small perturbations compared to the 
background permittivity e c , we obtain Q > and Q > P. 
In this way, vj > and v~ < 0. Summarizing, we obtain 
that = \v^\z and t^ - = — In other words, when 

the e-ray propagates in the z direction it will present a 
phase velocity which is generally different from the veloc- 
ity of itself propagating in the opposite direction (— z). 
Thus, for each direction q there will be only one value 
for the velocity of the extraordinary ray, which propa- 
gates anisotropically. Additionally, since the ordinary 
ray propagates isotropically, there will be birefringence 
phenomena in any direction but those where the velocity 
of both rays coincide. 

With the usual definition of the refraction index (n = 
1/v) we obtain the following result for the maximum 



difference between the propagation of both rays in the 
present situation 



nil — nj_ 



£3 



-— E 2 + —Je^JT c EB sin( 

£ c £ c 



~ {EB - y/e7JT c E 2 sine) cost 

Z8 C 



.(66) 



It is interesting to note that, by considering the wave 
vector in the opposite direction (—q), the result presented 
in Eq. (f66|) is modified as 



nil — nj_ 



-— E 2 - —JeHTcEBsmt 

£ c £c 



^- (EB + ^]T C E 2 sin 6) cos ( 



(67) 



In this case [Eq. ((67]) ] the terms in ei and £2 present the 
same sign. Thus, if £i/£2 > and 9 = ir/2 in the latter 
configuration of fields, the magnetoelectric contribution 
to the birefringence would magnify the effect. 

a. Crossed fields: The case of crossed electric and 
magnetic fields (9 = tt/2), which has been experimentally 
measured, is obtained directly from Eq. (f6"6"|) as 



II 



£ 1 ^2 



to 



-E z + — ^T C EB) . (68) 

£f! £n 



The term corresponding to the Cotton-Mouton effect 
does not appear since we have considered the particu- 
lar case of non- magnetic medium (fx = fj, c = constant). 
The inclusion of this term can be done simply by setting 
the permeability as a function of the magnetic field. 

A kind of magnetoelectric birefringence will still occur 
when the medium behaves as e = e(B) and both external 
fields E and B are present. In this case Eq. (|68p reduces 
to nil — nj_ — ii c £<2,EB, which shows that the magnetoelec- 
tric birefringence would be the dominant effect. In fact, 
it appears as the unique effect. This non-standard mag- 
netoelectric birefringence occurs also in the more general 
situation of non-crossed fields. In this case, the difference 
between the refraction indexes is given by 



nn — nj_ ~ fi c £2EB sin ( 



(69) 



We notice that its maximum value occurs when the 
fields are crossed (9 — tt/2). 

b. Jones birefringence: The Jones effect is another 
kind of magnetoelectric birefringence, also recently ob- 
served [2||, which occurs with parallel electric E and 
magnetic B fields perpendicular to the wave vector q. 
This case can be obtained directly from Eq. (|66p [or 
from Eq. (f6"T|) ] and results in 



-E 2 



£3 
2e c 



EB 



(70) 
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Magnetoelectric birefringence with crossed fields and 
Jones birefringence are very similar effects. Neverthe- 
less, they are independent to each other. Note that the 
terms in the product EB are proportional to different co- 
efficients. In the case of crossed fields it is proportional to 
the coefficient £2, which is coupled with the dependence 
on B 2 in the permittivity. On the other hand, in the case 
of Jones birefringence, it is proportional to the coefficient 
£3, coupled with the dependence on the product EB in 
the permittivity. 

Another interesting aspect related to these effects is 
the form of the front wave at a given instant of time. 
As we have shown, for all cases the e-ray propagates 
anisotropically. Particularly, the magnitude of the ef- 
fect is generally different in opposite directions. In order 
to visualize the effect of this anisotropy in the propa- 
gation direction, let us analyze a particular case where 
(EB) = 0, (qB) — and (qE) = cos-0. Let us con- 
sider a Cartesian coordinate system such that E = x 
and B = y, and solve the equation for the curves repre- 
senting the intersection of the normal surfaces with the 
AZ-plane. Since the o-ray propagates isotropically, its 
associated normal surface is represented by a circle in 
this plane, at any instant of time. On the other hand, 
the e-ray does not propagates isotropically and its nor- 
mal surface will be not symmetric with respect to the 
plane XY, as shown in Fig. [T] We notice the possibility 
to the velocity of the e-ray be greater or smaller than 
the velocity of the o-ray, depending on the direction of 
propagation. Further, depending on the parameters in 
the permittivity expansion, there will be more than two 
directions in which the velocities of both rays coincide. 

For any of the cases considered in this section the corre- 
sponding polarization vectors e M can be directly obtained 
by returning each specific solution in the set of equations 
(|22 )) - (|24)) and solving it for the coefficients a, b and c. 



V. CONCLUSION 

A tensorial formalism was here presented for the de- 
scription of monochromatic electromagnetic waves inside 
material media with nonlinear dielectric properties. The 
limit of geometrical optics was considered. The eigen- 
vector problem for general media was presented and 
solved for isotropic material media characterized by di- 
electric coefficients as e fJlU (E, B) and ^^(E, B). The so- 
lution for the wave propagation in isotropic media de- 
scribed by these coefficients was obtained directly from 
the general dispersion relation expressed by Eq. (|34|) . 
Using this main result, a detailed analysis of the bire- 
fringence phenomena was performed. Some well known 
effects were recovered and some new information was 
presented. The results include the description of the 
Kerr electro-optical effect, the Cotton-Mouton magneto- 
optical effect, and the recently measured magnetoelectric 
birefringences. 

It was shown that the Kerr birefringence in isotropic 




FIG. 1: Normal surfaces for the ordinary (circular solid line) 
and extraordinary (dashed line) rays propagating in an isotropic 
material media with dielectric coefficients given by e = e(E,B) 
and /j, = He — constant. The symbols e and o stand for the front 
waves associated with the e-ray and o-ray, respectively. The plot 
is based on Eqs. (|6ip and l|62p . where we have considered the 
field configuration as E = x and B — y, and also q-B = 0. The 
numerical values were chosen as to satisfy the approximations 
assumed in Section HV Bl 



media is fundamentally related to the dependence of the 
permittivity on the electric field as e(E). The applica- 
tion of an external magnetic field does not bring changes 
to the effect. Symmetrically, the necessary condition for 
the appearance of the Cotton-Mouton effect is the de- 
pendence of the permeability on the magnetic field as 
fJ-(B). Similarly, an applied electric field does not alter 
the results. 

In what refers to the magnetoelectric birefringence, 
some few remarks are in order. First, from the general 
dispersion relation it was found three possible solutions 
for the phase velocities. One is related to the propagation 
of the ordinary ray and the other two are related to the 
propagation of the extraordinary ray. The two generally 
distinct values for the extraordinary ray are related to the 
distinct values of its velocity in opposite directions, as it 
was explicitly shown for the case where q is perpendicu- 
lar to the plane that contains the electric and magnetic 
fields. For the case of parallel fields it was found that the 
propagation is symmetric in opposite directions. In both 
cases the propagation occurs anisotropically. Another re- 
markable result is the birefringence effect occurring with 
the permittivity e(B) and constant permeability /i c , in 
the presence of external fields. For this case, if the exter- 
nal electric field is absent there will be no birefringence. 
Nevertheless, if both the electric and the magnetic fields 
are present a kind of magnetoelectric birefringence ap- 
pears as the unique effect. This result shows that it is 
possible to produce such birefringence effect without any 
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other accompanying standard birefringence, as the Kerr 
or Cotton-Mouton effects. By symmetry, it can be in- 
ferred that a similar effect occurs with n(E) and constant 
permittivity. Finally, it should be stressed the behav- 
ior of the front wave associated with the extraordinary 
ray. Depending on each medium in which the propaga- 
tion occurs, the extraordinary ray can present its velocity 
greater or smaller than the velocity of the ordinary ray 
depending on the direction of propagation. Further, the 
coincidence between the both rays can occur in several 
directions. 



Appendix 

In this appendix, the Kerr effect is derived by means 
of the Cayley-Hamilton method 31]. Using the same 
assumptions presented in Section [IV A II we obtain, from 
Eq. (0, 

Z a T = (e - h a T - e'E<*E T - -^q a q T . (71) 
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The solutions of the eigenvalue problem stated by Eq. 
([3]) is formally given by det|Z Q T | = 0, which can be 
evaluated as 

(Z 1 ) 3 -3Z 1 Z 2 + 2Z 3 = 0, (72) 
with Zi (i — 1,2,3) representing the traces: 



Zi 
Z 2 

z 3 
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' q 2 


fl c U) J \ 


HcU 2 


3e'q 2 (qE) 2 






^ 6 



2e'{qEf 

HcLO 2 



e' 2 E 4 



e' 2 E\ 



2e'(gE) 2 

fJ, c UJ 2 



? A 



(73) 
(74) 



(75) 



Now, introducing the results lj73]l -([7S |l in Eq. ([72]). and 
using the definitions introduced in Sections IIIII and IIV1 
we obtain 



(1 - [i c ev 2 ) [i c e{e + e'E 2 )v 2 -e- e'E 2 {q ■ Ef 



= 0. 
(76) 



As one can see, the above equation presents two solutions, 
which are the same solutions obtained in Section TlV A II 
for the ordinary and extraordinary rays, as stated by Eqs. 
([40]) and (|4Tj) . Now, the Kerr birefringence is obtained 
in the same lines as discussed in the Section HV A 11 
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